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I. INTRODUCTION 



There is growing evidence that in addition to its ex- 
traordinary electronic properties^, graphene might also 
be an exciting material for spintronics, a technology that 
would be based on the spin of electrons rather than on 
their charge. The impetus to study spin-related phe- 
nomena in graphene comes from two directions: i) the 
experiments of Tombros et al (Ref . |3) showed that it was 
possible to inject spin into mono- and few layers graphene 
and measure spin signals in a spin- valve setup, and ii) the 
recent observation^il of band splitting in graphene due to 
spin-orbit interaction (SOI). Although the intrinsic SOI^ 
is expected to be weak in graphene (not exceeding^^^ 
~ 50/ieV), the breaking of the inversion symmetry by 
an external electric field or by the presence of a sub- 
strate can result in a substantial externally induced SOI. 
In particular, Varykhalov et al (Ref. (3) reported a spin- 
orbit interaction induced band splitting of « 13 meV in 
a quasi-free-standing graphene on Ni(lll)/Au substrate. 
The spin-orbit coupling was identified as Rashba-type 
SOl^-^° (RSOI) and it was attributed to the high nu- 
clear charge of the gold atoms that were intercalated be- 
tween the Ni substrate and the graphene layer to break 
the strong carbon-nickel bonds and make the graphene 
layer quasi-free-standing. The fact that gold intercala- 
tion can decouple graphene from the nickel substrate was 
also supported by density functional calculations^^ and 
that it may induce sizeable Rashba-type SOI was indi- 
cated by the computations of Ref. H. Furthermore, gold 
intercalation was used to decouple graphene grown on 
Ru(OOOl) substrate^^ where a band-gap opening at the 
Dirac-point of the graphene band structure was observed 
as well. The appearance of the gap was ascribed to the 
breaking of the symmetry of the two carbon sublattices in 
graphene. A gap opening in the graphene band structure 
was also found when the strong nickel-graphene bonds 
were passivated by potassium intercalation's,. Besides 



metal surfaces (for a review see Ref. IT^ . intensive re- 
search effort, both theoreticali^"— and experimental (see 
e.g. Refs. 3 , HgI. [T8l - l20l . earlier developments are reviewed 
in Ref. [2l|), is directed towards studying graphene on 
SiC substrate. These experiments show therefore that 
substrates can induce SOI and/or open a band gap in 
monolayer graphene. 

Angle-resolved photoemission spectroscopy (ARPES) 
is an important experimental technique that provides 
direct information on the bulk and surface electronic 
band structure of solid state materials [see e.g. Refs. [2^ 
and [2^. ARPES has also become a major tool to 
study graphene on various substrate a'^'^'^^d^'^^i^"'^"''' '— 
By also measuring the spin polarization of the photoelec- 
trons (the so called spin-resolved ARPES or SARPES 
technique'^^) and using a sophisticated data analysis 
method^ one may observe band splittings smaller than 
the intrinsic linewidth of regular ARPES experiments, 
providing a powerful tool to measure spin resolved elec- 
tronic bandstructure. Indeed, SARPES measurements 
were used in the experiments of Refs. H and to investi- 
gate the spin dependent band splitting in graphene. 

Our work is motivated by the fact that, as mentioned 
above, substrates can induce RSOI and/or open band 
gap in monolayer graphene. Therefore the interplay of 
the two effects, i.e. the RSOI and sublattice asym- 
metry induced band gap opening may be important in 
some systems. We note that small band-gaps (by which 
we mean a few tens of meV) are not easily detected 
by ARPES because of the finite experimental resolution 
(10 — 50meV) and because of the finite intrinsic linewidth 
in the measurements. We theoretically demonstrate that 
broken carbon sublattice symmetry coupled with RSOI 
induces a finite out-of plane spin polarization in mono- 
layer graphene, therefore SARPES measurements could 
detect small band gaps even if conventional ARPES can 
not. We study the constant-energy angular maps and the 
spin-resolved momentum distribution curves (MDC) '^^i'^^ 
of photoelectrons as a function of initial state energy and 
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line-broadening for finite RSOI and sublattice asymme- 
try induced band-gaps. Our work is therefore comple- 
mentary to Ref. |33| in which a similar study was pub- 
lished for zero SOI and also to Ref. 'si' which focused on 
the effect of RSOI on photoclcctrons but the sublattice 
asymmetry was not considered and the dependence of the 
MDC-s on initial state energy and line broadening was 
not discussed in details. 

The rest of the paper is organized in the following way: 
first, in Section |lT] we show that if both RSOI and sub- 
lattice asymmetry are present then the quasiparticles in 
monolayer graphene acquire a non-zero out-of-plane spin 
polarization in a part of the Brillouin zone. We then show 
in Section IIIII how the spin polarization of quasiparticles 
(both in and out-of plane) is related to the spin polariza- 
tion of photoclcctrons. Using these results in Section ITVl 
we present a numerical study on the initial-state energy 
and intrinsic line broadening dependence of fixed-energy 
ARPES angular maps and spin-resolved MDC's and we 
point out the signatures of sublattice asymmetry. Fi- 
nally, in Section |V] we discuss the possible experimental 
relevance of our work and give a short summary. Some 
details of the calculations in Section IIIII can be found in 
Appendix [X] 



II. RSOI AND SUBLATTICE ASYMMETRY IN 
GRAPHENE MONOLAYER 

In a previous publication^ we showed that starting 
from the tight-binding Hamiltonian suggested in Ref |5| 
to describe RSOI in monolayer graphene one can arrive 
at the following Hamiltonian in the continuum limit at 
the K point of the Brillouin zone (BZ): 
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[The BZ of graphene with the high symmetry points F, K 
and K' is shown in Fig.[lja).] Huso in Eq. ^ is written 
in the basis {\A f), \B \),\A i), \B |)} {{A,B} denoting 
the two triangular sublattice of graphene's honeycomb 
lattice and {f, i} is the basis in spin Hilbert space). The 
parameters appearing in the Hamiltonian ([T]) are as fol- 
lows: vf = 37oao/(2?i), where oq is the bond length be- 
tween the carbon atoms, 70 is the hopping amplitude be- 
tween next neighbour carbon atoms, v\ = SXaao/i^h), 
where A_r gives the strength of the RSOI in the tight- 
binding model of Ref '5'. Furthermore, p± — i: ipy 
and PxTpy SLie momentum operators. The Hamiltonian 
([T]) differs from the Hamiltonian put forward in Ref. [l^ 
by the terms v\p±. Like the terms ±3iAi?, they appear 
because of the spin-orbit interaction and they lead to 
trigonal warping of the bands at low energies, i.e. close 
to the K point of the BZ. Note that for wavenumbers 
far from K point there is another kind of trigonal defor- 
mation of the bands, which is a lattice effect [see e.g. in 
Fig.Hfb)]. It turns out that one can understand- - all the 



salient features of the spin polarization at low energies 
already without taking into account the v\p± term be- 
cause it gives higher order corrections in the wave vector 
k = {kx,ky) [measured from the K point]. Therefore in 
our analytical calculations we use the following Hamilto- 
nian: 
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where the terms ±A/2 account for a possible breaking of 
the symmetry of the sublattices A and B. 
The eigenvalues of Hamiltonian ^ are: 



e^,{k) = |y^4t;|,fi2fc2 + A2 + ISXj, - iyl8N{k) , (3) 
where fc = Ikl and 



(4) 



The index fi = 1(— 1) corresponds to conductance 
(valance) bands, whereas v — 1 for the low energy bands 
touching at k = (for A = 0) and u ~ —1 for the spin 
split bands^iiSii^i [a schematic of the band structure is 
shown in Fig. [Hb)]. In the case of AB asymmetry, i.e. 
for A ^ a gap opens in the spectrum at the Dirac-point 
(k-0). 

The RSOI leads to a particular spin polarization of the 
bands^i^. The expectation value of the three compo- 
nents of the quasiparticle spin in an eigenstate [^'^'''(k)) 
of the Hamiltonian ([2]) can be calculated as: 



' {Sx.y.z) — Tr (^Q^' Sx.y.z^ 



(5) 



Here Q^'^^ik) = |^'^''"(k))(*^''^(k)| is a 4 x 4 projector 
and a convenient way to calculate these projectors can 
be found in Appendix \^ The operator Sx,y,z is given 
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where I2 is the 2x2 identity 



matrix acting in the pseudospin space and <Jx,y,z are Pauli 
matrices acting in the quasiparticles' spin space. The 
expectation values of the spin components are found to 
be (in units of |): 



'"{Sx) 



-2hvFky\fi 
3JV{k) 



(Sy) 



2hvFkxXi 
3J\f{k) 



and 



{Sz 



2AA(fc)|e^,(k)| 



(6) 



(7) 



The X and y components of the spin polarization are 
independent of the sublattice asymmetry [Ens. (|6t] and 
we obtain the same results as in Refs. d, 113, S, and [H 
i.e. the in-plane component of the spin shows rotational 
symmetry around the K point, it is perpendicular to the 
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FIG. 1. a) Schematic of the hexagonal Brillouin zone of 
graphene with the T point and vectors K, K' pointing to 
the corresponding corners of the BriUouin zone. The coor- 
dinate system that we use in the momentum space is also 
shown, b) Schematic of the energy bands near the K point 
of the BZ, as obtained from Eq. @ for A = (solid) and 
A 7^ (dashed). The energy splitting between the spin-split 
bands is 3A_r for A = 0. If A 7^ 0, a band gap of A opens at 
k = 0. The (/i, 1/) indeces corresponding to a given band are 
also indicated, c) solid [dashed] line: expectation value of the 
z component of the spin as a function of k = |kt (in units of 
the carbon-carbon bond length ao) in the upper valance band 
(/i = -1, 1/ = 1) calculated from Eq. (0 [Eq. (O]. The inter- 
play of sublattice asymmetry and RSOI leads to a finite {Sz}- 
The width of the peak, eis defined in the figure, is independent 
of the asymmetry parameter A [see text below Eq. ((8)l]. 



wavevector k and its magnitude depends on k, vanish- 
ing at k = 0. One can see from Eq. (O that compared 
to the case of equivalent sublattices (A = 0) where the 
spin has only in-plane components for all bands^'ii^i^ii^, 
the interplay of sublattice asymmetry and RSOI leads 
to finite z spin polarization of electrons in the vicinity 
of the K point. For the v — 1 bands exactly in the K 
point (k = 0) the spins are fully polarized and perpen- 
dicular to the graphene sheet: ^'^{Sz) = n = ±1, while 
for the = — 1 bands the spin z component points into 
the opposite direction as in the v = 1 bands and it can 



be significantly smaller: ^'~^{Sz) = — yuA/y^^ 
(We note that ^''^{Sz) is the expectation value of the 
spin z component averaged over a unit cell and not on 
individual carbon atoms within the unit cell, which was 
discussed in Ref.[37|-) At the K' point, the other (inequiv- 
alent) point of the graphene BZ where the valence and 
conductance bands touch for A = 0, the spin polarization 
is exactly opposite than at the K point, as required by 
the time-reversal symmetry. Expanding the right hand 
side of Eq. ([7]) assuming that hvpk ^ Xr we find for the 
^' = 1 bands: 



^,1 



{Sz 



2h^vlk^ 



(8) 



range where the spin z component is non-zero. Fig. [Ijc) 
shows the z polarization computed using Eq. ([7|) and its 
approximation from Eq. (jS]). Estimating the width of 
the peak by the k values where Eq. becomes zero 
we find Afcs = SV^Xr/IIvf, which is independent of the 
asymmetry parameter A. Taking vp ~ lO^m/s and e.g. 
« 20meV we find that Aks ~ 0.01 A^^. 



III. THEORETICAL DESCRIPTION OF 
SARPES FOR GRAPHENE 

In this and the next section we will analyze the ef- 
fects of sublattice asymmetry on the SARPES spectra 
performing both analytical and numerical calculations. 
As in most of the relevant graphene literature^i^ii^S, we 
assume that the emitted photoelectrons can be character- 
ized by a simple plane wave of momentum p, spin a and 
energy Ep^a = (however, see e.g. Ref. Isol for the lim- 
itations of this assumption) . The flux of photoelectrons 
emitted from an initial state of momentum h(K. + k), 
energy epi,(k) in band {fi, v) is found to be 



'^P||/fi-(K-|-k+G),0 <5(^ + - Ep 
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[Some details of the calculations leading to Eq. (jH]) and 
Eq. (fH)) below are given in Appendix El] Here Q^^''(k) 
is a 2 X 2 projector onto the photoelectron spinor: 
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where Q^'"^ was introduced after Eq. ([5]) and the unitary 
matrix U is given by 



u= ( r- 
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which one can use to give an estimate of the wavenumber 



Here G is an arbitrary reciprocal lattice vector and r 
is a vector pointing from lattice site B to site A in the 
unit cell of graphene. In the following we will always 
take G = 0, since we will concentrate on one BZ. The 
Kronecker delta in Eq. (jl4p expresses momentum conser- 
vation (p|| is the component of the momentum of photo- 
electrons parallel with the graphene surface) . Finally, the 
Dirac delta function in Eq. (jl4p ensures the energy con- 
servation {W being the work function of graphene). We 
do not address the question of dynamical processes that 
lead to energy broadening but use a phenomenological 
approach by introducing a Lorentzian S{e) — >■ ^^'^pa (see 
Ref. mi) in the figures of Section lTVl with the parameter F 
representing the value of the broadening. To keep the for- 
mulae uncluttered we suppress henceforth the Kronecker 
and Dirac delta functions expressing the momentum and 
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energy conservation, they should be understood to ap- 
pear on the right hand side of Eqs. p^ - dTSl) below. 

Using the explicit form of the quasiparticle spinors, 
calculations detailed in Appendix [K\ yield 
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As in previous theoretical worksS^i^i^^ (which how- 
ever did not consider either RSOS^i^ or sublattice 
asymmetry^i) we find a strongly anisotropic photoelec- 
tron intensity [see e.g. numerical results in Fig. ^a)] 
which originates from sublattice interference^ and there- 
fore it is present^^ even if Aij = 0. Such anisotropy was 
observed experimentally^iii2ii^i2^ too. Indeed, Eq. 
for large wave numbers {hvpk ^ \^r\) can be approxi- 
mated by 
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from where it is easy to see that the intensity is mini- 
mal in the region where ky ^ \k^\^\^\l (tiv p) (J^v p) 
and fiky > 0. (By introducing the parametrization 
{kx,ky) — A; (sin 6*, cos 0) one can see that ^''^I takes a 
similar form to the result of Ref. [13, though in our no- 
tation the indices yit, v have slightly different meaning.) 
Since the intensity of photoelectrons tends to vanish in 
this region, the authors of Refs. and HI called this re- 
gion a dark corridor. For k — 0, in contrast, the intensity 
is isotropic. [We note that as the recent experiment of 
Gierz et al showed [Ref. (jsol )]. the angular distribution of 
photoelectron intensity also depends on the energy and 
polarization of the incident light. Our calculations should 
be relevant for p-polarized incident light.] 

In terms of Q^'^ , the expectation value of an operator 
O which gives the result of a measurement on photoelec- 
trons coming from band (/x, v) is given by 

''■"(0)(p) =Tr (^dQ'^''^(k)^ I Tr (q^'^'O^^)^. (14) 

We make use of Eq. (|14l) to calculate the photoelectron 
spin-polarization vector {'Px,'Py,'Pz) — jr{{sx), (sy), (sz)) 
where the operator s^^y^z = <^x,y,z acts on the photoelec- 
tron spin. Using Eqs. p4)) and p2|) we find for the com- 
ponents of the polarization that 

v\r (|A^ -I- Ivphky {vphky - £^^(k))) - f AflA^(fc) 
Mik)e^,,ik) ~ vphkyiN'ik) - 
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FIG. 2. a) intensity distribution; b) and c) spin-polarization 
of photoelectrons emanating from the upper valance band 
(/x = — l,i/ = l) for momentums close to the K point. In 
b) vector plot of the in-plane component of the spin polar- 
ization is shown. In c) density plot of the z component of 
spin indicates the regions in the BZ where the out-of-plane 
polarization is finite. We used A = 40 meV, \r — 66 meV in 
these calculations. 



It is immediately clear from Eq. (I15cp that similarly 
to Bloch electrons, photoelectrons also acquire a finite z 
polarization due to the interplay of sublattice asymmetry 
and RSOI. The magnitude of '^''^Vz is largest at the Dirac 
point for the 1^=1 bands where it reaches the value of 
unity. For the v = —1 bands the photoelectron polariza- 
tion is smaller: t'-^Vzik = 0) = -A/^A^ -|- 36A^. In 
fact, as the density plot in Fig. [5fc) shows for the upper 
valence band, ~^'^7^z(k) is finite everywhere in the dark 
corridor and is very small outside it (the plots for other 
bands are similar and thus not shown). This suggests 
that in a constant energy SARPES measurement the eas- 
iest way to observe the finite z polarization is to use en- 
ergies close to the Dirac point, otherwise one would have 
to collect data from the dark corridor, which is difficult 
due to the low photoelectron intensity and spin-detector 
efficiency. 

Regarding the in-plane component of the photoelec- 
tron spin, Ref. [s^ has found that in the case of equiv- 
alent sublattices it exhibits a rather peculiar behavior, 
especially in and close to the dark corridor, where the 
photoelectron spin is rotated with respect to the quasi- 
particle spin. Moreover, Ref. [13 also showed that the 
in-plane spin polarization of photoelectrons is not zero 
in the K point even though the mean spin of Bloch elec- 
trons is zero there [see Eqs. (]§])]. We find from Eqs. (|15ap 



and (jl5bp that the breaking of the AB symmetry does 
not alter significantly this picture of the in-plane polar- 
ization, thus we will only briefly discuss it. An example 
of the photoelectron in-plane spin polarization is shown 
in Fig. HJb) for the upper valence band. One can clearly 
observe that the spins are rotated in the dark corridor 
(at ky < 0, fca; 0, see Fig. HJa) where the intensity map 
is shown for the same band). In contrast to the in-plane 
spin of quasiparticles, the corresponding spin component 
of photoelectrons therefore does not show rotational sym- 
metry around the K point. 

The opening of a small gap at the Dirac point due 
to the AB symmetry breaking effect of a substrate is not 
easy to detect in an ARPES measurements because of the 
finite energy resolution of the experiments and because of 
the energy broadening of the bands. In the next section 
we investigate the possibility of detecting the sublattice 
asymmetry through photoelectron spin polarization. To 
this end we compute the intensity maps and spin polar- 
ization distributions of photoelectrons at given energies. 



IV. NUMERICAL (S)ARPES CALCULATIONS 

In this section we discuss the results of numerical cal- 
culations of constant-energy intensity maps and spin po- 
larizations along certain directions in the BZ. In Ref. 1 3 a 
we showed that the Hamiltonian of monolayer graphene 
for finite RSOI is formally the same as the Hamiltonian 
of bilayer graphene, if certain weak hopping amplitudes 
in the latter system can be neglected. The aim of this 
section is twofold. Firstly, we point out both the similar- 
ities and the differences in the constant energy ARPES 
intensity maps of monolayer graphene with RSOI and 
bilayer graphene. Secondly, we show photoelectron spin 
polarization calculations along certain directions in the 
BZ (spin-resolved MDCs) and relate them to the fixed 
energy ARPES intensity maps. In the calculation of spin- 
resolved MDCs we assume that the background is small 
and disregard its influence on the lineshapea^i. Since the 
AB asymmetry does not break the particle-hole symme- 
try of the Hamiltonian, we only show calculations for 
energies in the valence bands. We assume strong RSOI 
and use A^, = 66 meV corresponding to « 200 meV spin- 
splitting of the bands^. 

We start the discussion with intensity maps taken at 
energies close to the Dirac point. In the derivation of 
Eq. (IT^ we neglected those terms in the graphene Hamil- 
tonian which cause trigonal warping of the bands for 
low energies if RSOI is finite [see the discussion below 
Eq. ([IJ]. This approximation is useful to understand 
the main features of the spin-polarization but for strong 
RSOI the neglected terms do cause a noticeable change 
in the fixed-energy intensity maps. In the calculations 
shown below therefore we take these terms into account 
as well. 

In Figs, ^c) and ^d) only a small broadening of the 
lines is assumed. Because of the strong RSOI {Xn = 
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FIG. 3. Constant energy (S)ARPES calculations close to 
the Dirac point, a) and b): schematic band structure at the 
K point of the BZ for zero and finite sublattice asymmetry, 
respectively. Dashed lines indicate a constant energy cross- 
section a.t E = — 37meV where the photoelectron intensity 
maps in c) and d) are obtained, c) and d): constant energy 
intensity maps for A = and A = 40meV, respectively. The 
dashed lines indicate the direction in k space along which the 
spin polarization curves in e) and f) are calculated, e) and f): 
the X (dashed), y (dashed-dotted) and z (sohd) component of 
the photoelectron spin polarization for cross sections shown 
in c) and d), respectively, (left axis). The dotted lines indi- 
cate the intensity profile along the same cross-section (right 
axis) . Subfigures in the left [right] column correspond to sub- 
lattice anisotropy parameter A = [A = 40 meV]. Other 
parameters of the figure are A = 66meV and F = 12.5 meV. 



66meV), small F and low energy {E = — 37meV) the 
photoelectrons come predominantly from the upper va- 
lence band. One can observe the following important fea- 
tures: similarly to monolayer graphene with zero RSOI 
(Refs. m and ^8*) there is a characteristic angular vari- 
ation in the intensity which is due to sublattice inter- 
ference and particularly in Fig. [3][c) one can see a low 
intensity region (the "dark corridor") around ~ 
and ky < 0. Nevertheless, as a consequence of spin- 
pseudospin entanglement^ at these low energies the in- 
tensity distribution is more isotropic in the case of fi- 
nite RSOI than it is for A^; = 0. Fig. |31[d) shows that 
the main effect of finite sublattice asymmetry on the in- 
tensity maps is that it reduces the intensity anisotropy 
clearly seen in Fig. [3][c). Comparing e.g. Fig. [3Kc) and 
Fig. lUc) one can also notice that in the former figure 
there is a slight trigonal distortion in the intensity con- 
tour. This distortion, which is caused by the terms vxp± 
in the Hamiltonian ([!]), can only be seen for strong RSOI 
and close to the charge neutrality point. Note, that it is 
different from the trigonal distortion observable for en- 
ergies far from the Dirac point (see Fig. ^ which is a 
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lattice effect. 

Figs, ^e) and ^i) show tfie spin polarization as a 
function of momentum along the direction indicated by 
dashed line in Figs. [3jc) and[3l[d), respectively. As evi- 
denced by the peaks in the x polarization component Vx 
[and also noted in Ref.HJ], in contrast to Bloch electrons, 
the spin polarization of photoelectrons is not necesseraly 
transversal to the momentum k. One can also see that 
Vy changes sign as the fc^; = line is crossed. The out- 
of-plane component of the photoelectron spin is zero if 
no sublattice asymmetry is assumed [Fig. [3][f)] but Vz 
is finite if A 7^ 0, as in Fig. EJe). This means that 
through SARPES measurements in systems where RSOI 
is nonzero the AB asymmetry can be detected even if the 
sample is slightly p-doped, i.e. states around the Dirac- 
point are not directly available by ARPES. 

In Fig. |4] the constant energy maps are calculated at 
E — — llOmeV, i.e. not in the close vicinity of the charge 
neutrality point. As the schematic figures |ll[a) and ^h) 
show, because of the large spin-splitting (and a small 
broadening of F = 16.7 meV) assumed, all the photo- 
electrons would still originate from the same band as in 
the previous case. The intensity maps in Figs. Hl[c) and 
(d) resemble closely the corresponding maps of mono- 
layer graphene (see e.g. Fig. 2 in Ref. Issl ). In particu- 
lar, one can observe an almost complete suppression of 
intensity in the dark corridor and the disappearance of 
the trigonal distortion of the intensity maps, apparent 
in Figs. ISKc) and (d). Furthermore, comparing Fig. |31[c) 
and Fig. SKd) one can see that the presence of a small 
asymmetry gap [A = 40 meV in Fig.UKd)] would be prac- 
tically undetectable in an ARPES measurement at this 
energy. Nevertheless, as Fig. |31[f) and Fig. SKh) show, if 
A 7^ there is a small but finite z polarization compo- 
nent. Comparison of Fig. and Fig. |ll[h) illustrates 
the feature shown in Fig.[2][c): 7^2(k) is largest in the dark 
corridor, therefore in a constant energy measurement it 
is larger if the direction in the k space is chosen such 
that it is closer to the dark corridor [ Fig. is calcu- 
lated for ky = — O.OOTA"^ with maximal polarization of 
pmax ^ _o.i2, whereas ky = -QAAk'^ in Fig.gKh) and 
pm.ax ^ _o.i7j. Note, that even in the case of Fig. EKh) 
the curve is not actually calculated in the dark corridor, 
the ARPES intensity peaks (shown by dotted line) for 
this cross-section are roughly 40% of the maximum in- 
tensity that can be found at this energy [black arc close 
to the upper edge of Fig. lU^d)]. On the other hand, the 
in-plane components of the spin polarization are practi- 
cally the same for the A = [Figs. IDJe), IDJg)] and A 7^ 
[Figs.|4i;f),|4i;h)] cases. 

In ARPES measurements the energy broadening is of- 
ten quite substantial. To see the effects of broadening on 
the SARPES spectra we repeated the calculations shown 
in Fig. 2] for a larger broadening parameter. The results 
for F = 83.5 meV are presented in Fig. [5j Although the 
ARPES fixed energy contours are significantly blurred 
due to the large F, [Figs, ^c) and [SJd)] the broaden- 
ing would actually lead to a bigger out-of-plane spin po- 




FIG. 4. Constant energy (S)ARPES calculations for low ener- 
gies and small broadening F. a) and b): schematic band struc- 
ture at the K point of the BZ for zero and finite sublattice 
asymmetry, respectively. Dashed lines indicate a constant en- 
ergy cross-section at _E = — 110 meV where the photoelectron 
intensity maps in c) and d) are obtained, c) and d): constant 
energy intensity maps for A = and A = 40 meV, respec- 
tively. The dashed line at ky = -O.OOTA"^ [ky = -0.014A"^] 
indicate the direction in k space along which the spin polar- 
ization curves in subfigures e) and f) [ g) and h) ] are cal- 
culated, e) and g) [ f) and h) ]: the x (dashed), y (dashed- 
dotted) and z (solid) component of the photoelectron spin 
polarization for the two cross sections shown in c) [ d)] (left 
axis). The dotted lines indicate the intensity profile along 
the same cross-section (right axis). Subfigures in the left 
[right] column correspond to sublattice anisotropy parame- 
ter A = [A = 40 meV]. Other parameters of the figure are 
A = 66 meV and V = 16.7 meV. 



larization amplitude see Figs. (S^f) and Efh)], hence it 
would make the detection of the z spin polarization eas- 
ier [c.f Figs. Hff) andSlJh)]. This happens because for 
large broadening electrons having energies closer to the 
Dirac point can also contribute and they have larger spin 
z component. Other noticeable feature in Figs.[ni[e)-[SJh) 
compared to Figs. IHe)-|ll[h) is that one can clearly see 
that Vy changes sign three times for small kx values. This 
is not apparent in e.g. Figs.HJe) and|4l[f) because of the 
small amplitude of these oscillations there. 

Finally, we consider the constant-energy intensity 
maps and spin polarizations at energy E = —660 meV, 
i.e. quite far from the Dirac point. For these calcula- 
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FIG. 5. The effects of the broadening parameter T on the 
(S)ARPES spectra. As in Fig. [4l all calculations are for 
E = — llOmeV. a) and b): the same as in Fig. O c) and 
d): constant energy intensity maps for A = and A = 
40 meV, respectively. The dashed line at ky — — 0.007A~^ 
[ky = — 0.014A~^] indicate the direction in k space along 
which the spin polarization curves in subfigures e) and f) [g) 
and h)] are calculated, e) and g) [ f) and g) ]: the x (dashed), 
y (dashed-dotted) and z (solid) component of the photoelec- 
tron spin polarization for the two cross-sections shown in c) [ 
d) ] (left axis). The dotted lines indicate the intensity profile 
along the same cross-section (right axis). Subfigures in the left 
[right] column correspond to sublattice anisotropy parameter 
A = [A = 40 meV] . Other physical parameters of the figure 
are A = 66 meV and F — 83.5 meV. 



tions we used the tight-binding Hamiltonian of Ref. [H. 
Since this energy is larger than the spin splitting SXr = 
200 meV used in our calculations, both valence bands 
contribute to the ARPES and SARPES spectra. We as- 
sume for simplicity that the broadening F is the same 
for both bands and present calculations with two differ- 
ent Fs, the first one being much smaller than the spin- 
splitting of the bands, while the second one is comparable 
to it. As » A in this case, the ARPES and SARPES 
spectra are practically the same for A = or A 7^ 0, 
therefore we only show results for A = 0. 

If the broadening is moderate, as in Fig. ^h), there 
are two discernible ringlike patterns, each corresponding 
to photoemission from states in one of the two bands. 
The rings show slight trigonal distortion, but in contrast 
to Fig. [2Uc), this is a lattice effect and would also be 



FIG. 6. Constant energy (S)ARPES calculations far from 
the Dirac point, a): schematic band structure at the K point. 
Dashed line indicates a constant energy cross-section at i5 = 
—660 meV where the photoelectron intensity maps in b) and 

c) are obtained, b) and c): constant energy intensity maps for 
F = 50 meV and F — 134 meV, respectively. The dashed lines 
at ky = 0.035A~^ [ky = — 0.035A~^] indicate the direction in 
k space along which the spin polarization curves in subfigure 

d) and e) [ f) and g) ] are calculated, d) and f) [e) and g)]: 
the X (dashed), y (dashed-dotted) and 2 (solid) component of 
the photoelectron spin polarization for the two cross sections 
shown in b) [ c) ] (left axis). The dotted lines indicate the 
intensity profile along the same cross-section (right axis). 



observable^ for = 0. The double ringlike pattern 
is reminiscent of the intensity maps found for bilayer 
graphene at high energies^^*^, but an important differ- 
ence is that in Fig. both rings have approximately 
the same intensity. The similarities between the ARPES 
maps of the two systems are due to the similar bandstruc- 
tures (for a discussion of the relation between the Hamil- 
tonians of monolayer graphene with RSOI and bilayer 
graphene see Ref. |36j). The difference in the intensity 
patterns stems from the fact that there are four carbon 
atoms in the unit cell of bilayer while there are only two 
in monolayer graphene therefore the transition matrix el- 
ements in the photoemission calculations are different. 

If the broadening is substantial, as in Fig.lHl^c) , the two 
rings are no longer easily discernible (and they may even 
completely overlap). Nevertheless, as the dashed-dotted 
curves in Figs.lHlJe) andlHIg) demonstrate, the y compo- 
nent of the spin polarization changes sign as a function 
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of kx roughly in the middle of the intensity peak (dotted 
line). This is an indication that two bands are involved 
in the photoemission, as the sign of '^''"Py is different 
for the v = \ and v = —1 bands [see Eq. (|15bp ]. Fur- 
thermore, comparison of Figs. [6je) and [Hg) [Figs. El^d) 
and [6jf )] shows that the overall shape and the number 
of sign changes in ^'^T^y do not depend on whether it 
is calculated for a positive or negative ky = const value 
[see Fig. [Ifb) or ^c) for the cross-sections along which 
Figs-[nid)-[njg) were obtained]. In contrast, for ^'"Vx the 
number of sign changes in the low intensity region (small 
|fca;| values) is affected by the choice of the ky, as e.g the 
comparison of Figs. |6l[e) and|6l[g) can illustrate. 



V. DISCUSSION AND SUMMARY 

We would first briefly comment on the experimental 
relevance of our results. As mentioned in the Intro- 
duction, a significant spin-orbit coupling was found in 
gold intercalated Ni(lll)/graphene system^ and the SOI 
was attributed to the presence of the gold atoms. Spin- 
resolved MDCs were not shown however in Ref.j^. Subse- 
quently, Ref . vl proved that gold intercalation can decou- 
ple graphene from the Ru(OOOl) surface as well. Another 
notable recent development is that gold intercalation has 
also been used for the Si face of SiC substrateAS' where 
the strong covalent bonds between the SiC (0001) and 
the first graphitic layer were suppressed by this method 
resulting in a slightly p-doped graphene that was only 
weakly infiuenced by the substrate. SARPES measure- 
ments were not published in Refs. [l^ and [l^, though 
it would be interesting to know if gold can induce SOI 
in these systems as well. The Ru(0001)/gold/graphene 
system appears to be particularly interesting from our 
point of view because ARPES measurement indicate a 
band gap > 100 meV, so that if RSOI is non-zero in this 
system then a finite out-of-plane polarization of photo- 
electrons should be measurable. A qualitatively similar 
polarization to the one predicted by this model, with an 
"abrupt rotation of the spin" at the K point of the BZ 
was measured when thallium was deposited on Si(lll) 
surface^, though Ref. |4l] explained the effect by the pres- 
ence of a local effective magnetic field. Finally, we note 
that Ref. reported a large and anisotropic spin splitting 
in graphene, including a nonzero out-of plane polariza- 
tion component, but the origin of this effect seems to be 
unclear at the moment. 

In summary, we studied the effect of RSOI and sub- 
strate induced sublattice asymmetry on the spin polariza- 
tion of quasiparticles and of photoelectrons in graphene. 
The breaking of AB sublattice symmetry opens a gap in 
the band structure of graphene at the K point of the BZ. 
If RSOI is finite, the interplay of the two effects induces 
a non-zero out-of-plane component of spin polarization 
of quasiparticles in part of the BZ. RSOI also affects the 
intensity and spin distribution of photoelectrons, hence 
it can be studied with the (S)ARPES technique. For 



strong RSOI, the fixed-energy intensity maps taken at 
low energies, close to the K point of the BZ, show a char- 
acteristic trigonal deformation. This deformation of the 
intensity map survives the switching-on of an AB symme- 
try breaking potential given by the asymmetry parameter 
A, as long as A is much smaller than the RSOI induced 
band splitting. Our spin-resolved MDCs calculations also 
show that an important sign of the simultaneous presence 
of RSOI and sublattice asymmetry is if non-zero out-of- 
plane photoelectron spin polarization can be measured. 
It is important however, especially if A and RSOI are 
small, to choose the energy at which the spin-resolved 
MDCs are taken as close as possible to the Dirac-point, 
because for energies far from it the out-of-plane polar- 
ization remains finite only in the " dark corridor" , where 
the low photoelectron intensity would hinder the obser- 
vation of this effect. A carefully chosen cross-section in 
the momentum space or a large intrinsic energy broad- 
ening may, however, facilitate the observation of the spin 
z polarization in MDCs even at higher energies. Mean- 
while, the in-plane components of photoelectron polariza- 
tion remain qualitatively the same regardless of whether 
A is zero or not. If the fixed-energy intensity map is ob- 
tained at energies larger than the energy separation of 
two spin-split bands and their intrinsic energy broaden- 
ing r is small compared to their RSOI induced energy 
splitting, then the resulting ARPES calculation shows a 
double ring-like structure. For large F, the two rings may 
not be discernible any more, but SARPES measurements 
can nevertheless reveal the true band structure because 
of the sign-changes in the polarization components. 
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Appendix A: Outline of the theoretical SARPES 
calculations 

Here we briefly describe the calculation leading to 
Eq. (|14p . The Hamiltonian of the interaction between the 
Bloch electrons and the electromagnetic field in dipole 
approximation'^^ is given by 

Hir^t oc -^AV , (Al) 
1 

where A = Aqc'^'''""'^*^ is the vector potential. The tran- 
sition probability between an initial Bloch electron state 
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k, i^)) and a photoelectron state |p,<t) will be pro- p|| is the projection of momentum p onto the plane of 
^r,r^-^^,T^al fr, I /' F7. .\P'°' ^( h,., J- £ ^^(l^^ — Ep ,j — W) , graphene. By defining 



portional to | (i?jnt)k,'Qx,j.) ? ^(^^ 
where the photoexcitation matrix element is 



^{p,a\H,nt\K{^^,l^))■ (A2) 



p,cr; 



(A5) 



Explicit expression for (^^mt)k(^jy) can obtained by as- 
suming that the wavefunction of a photoelectron given 
by a plane wave |p,cr) oc e'P''/''|cr) and the wavefunction 
of a Bloch electron is 



the expectation value of an operator O with respect to 
the photoelectron state emanating from an initial Bloch 
state of momentum h(K. + k), energy ep^(k) in band 
(/i, v) can be calculated as 



|k, = 



1 



(A3) 



Here are vectors pointing to sublattice sites j — 
{A, B} in unit cell n, N is the number of unit cells in the 
sample, ipj^, (k) are the amplitudes of Bloch electrons on 
sublattice j with momentum k and spin a' and finally, 
<i>(r) is a Pz atomic orbital. The photoexcitation matrix 
element then reads: 



7V$p(Ap)U^;,(k)-He'°-V^^,(k) 



(A4) 

In Eq. (|A4[) $p is the Fourier transform of the atomic or- 
bital $(r), and G = mibi -I- m2b2 is a reciprocal lattice 
vector which is given in terms of primitive lattice vec- 
tors bi — {2n/ao,2Tr/^/3ao), b2 = (27r/ao, — 27r/\/3ao) 
and integers mi, m2. Furthermore, t = — R^ and 



^Pii/;i-(K+k+G).o ^(fii^ + e^i,{k) - Ep,^ - W). 

(A6) 



-(0>(p) 



Eq. (dH) then follows from Eqs. and (|X6|) . We 

note that a convenient way of calculating the projectors 
(5^'''(k) which are necessary to evaluate Eqs. ^ and (fT4)) 
[see Eq. pH]) ] is to make use of the following: if one de- 
notes hy Ej , j — 1 ... rid, Ud < N the distinct eigenvalues 
oi a N X N hermitian matrix H, then the projector onto 
the ryth eigenstate is given by the expression 



n( 



H-E,I 



(A7) 



n (Er, 

which does not necessitates the calculation of the eigen- 
vectors. In the mathematical literature the projectors 
are known as Frobenius covariants*^. In terms of the 
projectors and eigenvalues Ej the matrix H is given 
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